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The structural and electronic properties of germanium 
clathrates Ge46 and KsGe46 are studied by first principles 
calculations within the local density approximation. The 
equilibrium structures are obtained by ab initio pseudopoten- 
tial calculation combined with dynamic minimizations. The 
clathrate structure is found as a low energy phase for germa- 
nium. The electronic band structures for Ge46 clathrates are 
calculated and the band gap is found to be considerably larger 
than that of the diamond phase. Due to the effect of pentag- 
onal rings, strong similarity in electronic properties between 
clathrate and Ge24 fullerene structure are found. The effect of 
doping clathrate cages with metal atoms are examined. The 
KsGe46 clathrate is found to be metallic with the conduction 
bands only slightly modified by K dopants. 

71.20.Tx, 61.48.+C 



I. INTRODUCTION 

The Si and Ge clathrates can be viewed as covalent 
fullerene solids which is composed of three-dimensional 
networks of fullerene cages connected by face sharing. 
The silicon clathrate compounds Ma;Si46 and M^Siiac 
(M=Na, K, Rb, and Cs) were first synthesized in 1965E1. 
The structures of semiconductor clathrates can be clas- 
sified into two cubic structural types with 46 atoms or 
136 atoms per unit celllil. As shown in Fig.l, the type I 
clathrate (Clathrate-46) is formed by two smaller 14-face 
pentagonal dodecahedra (12 five-fold rings, Ih) and six 
larger 16-face tetrakaidecahedra (12 five- fold rings and 
2 six-fold rings, D 6( j). The structure of type II clathrate 
(Clathrate-136) consists of sixteen smaller pentagonal do- 
decahedra (12 five-fold rings, Ih) and eight larger hex- 
akaidecahedra (12 five-fold rings and 4 six-fold rings, T^j. 
Electronic structure calculations based on different meth- 
ods have shown that these two ocen network structures 
have similar electronic properticsou. 




FIG. 1. Atomic structure of fully-relaxed KgGe4g clathrate. 
Lattice constant of the simple cubic cell is 10.45 A. 

The research interests of semiconductor clathrates 
come from several aspects: (1) possible alteration 
of electronic structures and energy gap from stan- 
dard diamond formotj, (2) metal-insulator transition 
in M^Suae with different concentration of metallic 
impurityoij, (3) the, finding of superconductivity behav- 
ior in Na a; gaj ) Si46lj, (4) candidate for thermoelectric 
applicationsQ, (5) template to from three dimensional 
arrays of nanosized clustersu, (6) similarity in struc- 
tural and electronic properties between semiconductor 
clathrates and nanoclustersBEl 

There are lots of the experimental and theoretical Stud- 
ies on Si46 and Sii36 clathrates and their compoundsErO. 
The structural, electronic and vibrational properties of 
Si clathrates are investigatedpaLiJarious theoretical kfe 
els ranging from ab iniiica'0OE3 to tight-binding3'E3 
and empirical potentials. The most interesting result 
from those calculations is that the band gap is about 
0.7 eV higher than that of diamond phase. Experi- 
mental works |-,on Si clathrates include resistivity and 
magnetizations, transport propertujal, photoemission 
spectroscorayEI, NMRMrtll, RamanBH, ESRE3, neutron 
scatteringnil etc. 

In contrast to the intensive studies on silicon 
clathrates, our current knowledge on germanium 
clathrates is rather limited. Recently, the germa- 
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nium clathrate compounds such as K 8 Ge46, Rb a: Ge46, 
Naa;Gei36, CsgNai6 have been synthesized and-, their 
structures are analyzed with X-ray diffractioncatl An 
empirical potential calculation has been performed on 
pure germanium clathratestJ. However, there is no first 
principles electronic structure calculation and their fun- 
damental electronic properties are still unclear theoret- 
ically. In this work, we report results of first principles 
study on the structures and electronic properties of Ge46 
and KgGe46 clathrates. The equilibrium structure, elec- 
tronic band and band gap, electronic density of states and 
electron density distribution are obtained and discussed. 



type I clathrate is 15.3 % and its energy is 0.071 eV per 
atom higher than diamoncO. 

Table I The structural properties and band gaps of per- 
fect and relaxed clathrate compared with previous LDA 
local orbital (LDA-LO) and tight-binding molecular dy- 
namics (TBMD) calculation on Si46 clathrate. V/Vo is 
the reduced atomic volume. I denotes the Ge— Ge bond 
length, a is the equilibrium lattice constant. E g is the 
electronic band gap (indirect gap for both diamond and 
clathrates^. The numbers in brackets are experimental 
valuea23'E3. 



II. COMPUTATIONAL METHODS 

First principles SCF pseudopotential method&S is 
used to perform static calculation on the electronic 
structures and total energy of Ge46 and K 8 Ge4 6 in 
ideal clathrate structures with different lattice con- 
stants. The ion-electron interaction is modeled by numet 
ical BHS norm-conserving nonlopal pseudopotentialEZI 
in the Klcinman-Bylander formlHa. The Ceperley- 
Alder's exchange-correlation parameterized by JSerdew 
and Zunger is used for the LDA in our programed. The 
kinetic energy cutoff for plane-wave basis is chosen as 12 
Ryd. Ten symmetric k points generated in Monkhorst- 
Pack spiritEj are employed to sample the Brillouin zone. 

From the static calculation, the equation of states with 
ideal crystal structures and the equilibrium lattice con- 
stants are obtained. These ideal structures are further 
optimized by using structuraloainimization via conjugate 
gradient technique (CASTEPEl). The CASTEP program 
is based on plane-wave pseudopotential technique. It can 
relax the atomic position by computing the force acting 
on atoms from electronic calculation and moving atoms 
efficiently in a numerical wayc3. 





V/Vo 


l(k) 


a (A) 


E g (eV) 


Ge Diamond 


1.00 


2.407 


5.56 (5.66) 


0.40 


Ge46 (ideal) 


1.148 


2.375, 2.540 


10.43 


1.25 


Ge46 (relaxed) 


1.128 


2.380 ~ 2.433 


10.37 


1.46 


KgGe46 (ideal) 


1.271 


2.562, 2.739 


10.79 




KgGe46 (relaxe^) 


1.155 


2.346 ~ 2.492 


10.45 (10.66) 




Si 46 (LDA-LGOo 


1.17 


2.38 - 2.43 


10.35 


2.50 


Si 46 (TBMD)H 


1.14 


2.367 


10.20 


1.82 



III. THE STRUCTURES AND BAND GAPS OF 

GE46 



In Fig. 2, we present the equation of states for both 
Ge46 clathrate and diamond phase that is obtained from 
first principles static calculations. We find that the Ge46 
clathrate is a locally stable structure and its energy is 
only about 0.08 eV per atom higher than that of dia- 
mond phase. For comparison, the /?-tin phase of germa- 
nium is, about 0.25 eV higher in energy than the diamond 
phaseEx The low energy feature of clathrate phase is a 
natural consequence of its four coordinate characteristics 
and may be partially attributed to the softness of the 
bond-bending distortion modesa. As compared to dia- 
mond structure, the volume per atom in the clathrate 
phase is increased by about 14.8%. All of these results 
are very close to the previous empirical potential simu- 
lation on Ge46 , in which the change of atomic volume in 
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In previous theoretical study on bulk silicon and ger- 
manium of various phasesE_l, considerable similarity in 
the bonding behavior and phase diagram are found be- 
tween silicon and germanium. In Table I, we summarize 
our results on the structural properties and band gaps 
for perfect and relaxed Ge46 and KsGe46 clathrate and 
compare them with the previoujLLDA and tight-binding 
calculations on Si46 clathratealim. We find the difference 
of volume and energy between clathrate and diamond 
phase for germanium and silicon are comparable. 
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FIG. 2. Energy per atom of diamond and Ge46 clathrate 
as a function of atomic volume. 

The ideal clathrate structure for the Ge46 at equilib- 
rium lattice constant is further optimized with CASTEP 
plane-wave pseudopotential calculation. Both atomic po- 
sitions and unit cell parameters are allowed to relax. The 
SCF pseudopotential code used in static calculation has 
also been used to test the total energy for the initial and 
final structure and the results agree with CASTEP calcu- 
lations. After optimization, the lattice constant of simple 
cubic unit cell decreases from 10.43 A to 10.37 A. The 
relative atomic positions are only slightly relaxed from 
their initial configuration. The range of bond length dis- 
tribution has also been narrowed, i.e., from the 2.375 A 
~ 2.540 A in ideal clathrate structure to 2.38 A ~ 2.433 

A. 
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FIG. 3. Electronic band structure for (a) Ge46, (b) KsGe46 
clathrate near Fermi level. 

The electronic band structure calculations are carried 
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out for Ge46 clathrate with the equilibrium structure. 
The near band-gap band structures for the minimum 
energy configurations of Ge46 clathrates are plotted in 
Fig. 3 (a). Both the valence band maximum and conduc- 
tion band minimum are located on the r to X line and 
they are very close in k space. An indirect gap of 1.46 eV 
is found. As for bulk germanium in diamond phase, our 
calculations have yielded an indirect band gap of 0.40 eV 
between L and T points. The underestimation of band 
gap is a common feature of LDA calculations. Neverthe- 
less, our current results suggests that the band gap of 
Ge46 phase is about 1 eV higher than that of the dia- 
mond phase, which is similar to the 0.7 eV incpenient in 
band gap found for Si46 and Si34 clathratesa#BlI9. Fur- 
thermore, the electronic band structure of Ge46 clathrate 
are calculated independently via a simple tight-binding 
model. The tight-binding hopping parameters for ger- 
manium are taken from Ref.[32] and the parameters for 
nearest neighboring atoms in clathrate structure is as- 
sumed to be the same as that in diamond. According 
to tight-binding calculation, the band gap of diamond is 
1.13 eV and the gap for relaxed clathrate Ge46 is 2.93 eV. 
Although the tight-binding method usually overestimate 
the band gap of a system, the increment of 1.80 eV in 
the band gap from tight-binding model is in reasonable 
agreement with ab initio calculation. Therefore, we ex- 
pect that the true increase of band gap A from diamond 
to Ge46 clathrate phase might be between the LDA and 
TB prediction, i.e., 1.06 eV < A < 1.80 eV. This result 
suggests that the clathrate materials might be useful in 
the new electronic and optical application in the future. 

Since the Ge46 clathrate structure is essentially a 3D 
network composed by Ge2o and Ge24 cages with face 
sharing, it is interesting to compare the electronic prop- 
erties of the individual Ge2o and Ge24 cages with that 
of the clathrate. SCF pseudopotential electronic struc- 
ture calculations on isolated Ge2o and Ge24 clusters in 
fullerene cages are performed. The clusters are placed 
in a large simple cubic supercell with length of 28 a.u. 
In Fig. 4, we present the calculated electronic density of 
states (DOS) for Ge2o and Ge24 cages along with DOS 
for diamond and Ge46- The detailed analysis are given 
in the following. 
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FIG. 4. Electronic density of states (DOS) of Ge2o, Ge24 
fullerene, bulk germanium in Ge46 clathrate and diamond 
phases Gaussian broadening of 0.136 eV is used. 

Firstly, we find that most of the peaks in the DOS 
of Ge46 can be assigned to Ge24, while the DOS of 
Ge2o also show some similarity with Ge46- The DOS of 
Ge24 is closer to that of clathrate because most atoms in 
clathrate are associated with Ge24 cage. The presence of 
s-p gap in due to the large number of five-member rings in 
the cage and clathrate structure. This remarkable simi- 
larity implies that the Ge46 clathrate can be taken as the 
small Ge fullerene assembled solid with both geometri- 
cal and electronic hierarchy. However, both the Ge24 
and Ge2o cages are not semiconductor systems and do 
not have open band gaps like Ge46 clathrates. The band 
splitting in Ge46 can be attributed to the sharing of face 
atoms by Ge2o and Ge24 and the interaction between 
neighboring fullerene cages. 

On the other hand, we can compare the DOS of 
clathrate and diamond in Fig. 4. Several significant dif- 
ferences can be found. Besides the improvement of 
band gap in clathrate upon diamond phase, the total 
width of valence band of Ge46 (about 12 eV) is narrower 
than that of diamond (about 13 eV). This phenomenon 
has been also predicted in Si46 clathrateO and observed 
experimentallytd. We also found a gap opening between 
the s-like and p-like states in the case of Ge46 clathrates 
as well as Ge2o, Ge24 clusters, while s and p states are 
overlapped in the DOS of diamond. These remarkable 
differences in the electronic structure of clathrate and 
diamond can be understood by the large portion of fivc- 
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fold ring in clathrate structure. In the diamond lattice 
consisting of 100% six-fold ring, the 4s orbital can form 
complete antibonding states, which can distribute to high 
energy and overlap with the lower 4p-like states. In con- 
trast to diamond, the Ge46 clathrate is composed of 87% 
five-fold rings and 13% six-fold rings. In the five-fold ring, 
4s orbital cannot form complete antibonding states so 
that the top of 4s-like states is still lower in energy than 
the bottom of 4p-like and a gap between s and p states 
inside valence band opens. The similar effect of five-fold 
ring on the Ap bonding orbitals will induce the incom- 
pleteness of 4p-like states. As a consequence, the valence 
band top of Ge46 is lower than that of diamond, which 
corresponds to the narrowing of valence band width and 
broadening of fundamental gap between valence and con- 
duction band. The novel electronic properties caused by 
five-fold ring and the similarity between clathrate and 
fullerene cages may be explored for future electronic and 
optical applications. 

IV. ELECTRONIC STRUCTURES OF K 8 GE 46 
CLATHRATE 

Although the pure silicon and germanium semiconduc- 
tor clathrates are predicted to be locally stable, exper- 
iments have only synthesized the metal-doped clathrate 
compounds of Si and Ge. The metallic impurity atoms 
inside the clathrate might influence the electronic proper- 
ties of clathrate due to the interaction between the metal 
atoms and semiconductor skeleton. Here we have chosen 
the KgGe46 clathrate as such a model system to study 
the doping effect on germanium clathrate. 
The structural properties of K 8 Ge46 ar e studied with the 
same procedure we have applied on Ge46 clathrate: a 
static calculation of equation of states followed by a full 
relaxation. The minimized structure of K 8 Ge46 clathrate 
is shown in Fig.l and the structural parameters of this 
minimum energy configuration is given in Table I. It is 
natural to find that the lattice constant of cubic unit cell 
for K§Ge46 clathrate is larger than pure Ge46 clathrate in 
both ideal and relaxed structure because of the K atoms 
inside the fullerene cage. We also find a decrease of lat- 
tice constant from 10.79 A to 10.45 A after structure 
minimization. The minimized lattice constant (10.45 A) 
for the K§Ge46ris in good agreement with experimental 
value (10.66 A)E3. During the minimization, the geome- 
tries of KgGe46 relax from the perfect clathrate network 
a little more than that happened in the Ge46- 

First principles SCF pseudopotential electronic struc- 
ture calculation has been performed on the relaxed 
KgGe46 clathrate. In Fig. 3(b), we presented the elec- 
tronic band structure of KsGe46 near Fermi level. The 
system is found as metallic due to the K dopants. We 
can further examine the highest valence bands and lowest 
conduction bands of Ge46 and K§Ge46 shown in Fig. 3(a) 
and (b) in detail. The valence bands of KsGe46 are very 



closed to the original bands in Ge46, while the conduc- 
tion band structures have been slightly modified upon 
the inclusion of K atoms. 




Energy (eV) 



FIG. 5. Electronic density of states (DOS) of Ge46 and 
KgGe46 clathrate. Gaussian broadening of 0.136 eV is used. 
Note the similarity in the two DOS. 

We also study the difference of the electronic proper- 
ties between the pure and doped systems by comparing 
their the electronic density of states (DOS). The density 
of states for Ge46 and K§Ge46 clathrate are compared 
in Fig. 5. The DOS of valence electrons in KsGe46 is 
very close to that in Ge46 while the DOS for conduction 
electrons of KgGe46 ar e somewhat different from that 
in Ge46. On the other hand, the gap between valence 
and conduction band is 0.23 eV narrower in the case of 
K 8 Ge 46 . 
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complete charge transfer from K sites to Ge frames is 
found in the KsGe46 clathrate. 




FIG. 6. Contour plots of the electron densities of KgGe46 
on the (100) plane. No charge density is found at K site 
(marked by X), indicating complete charge transfer. 
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The analysis on the charge transfer and chemical bond- 
ing effects have been presented in the contour plots of 
the electron densities of K 8 Ge46 on the (100) plane in 
Fig. 6. We have also calculated the charge density distri- 
bution of Ge46 and find it is very close to that of K 8 Ge4 6 . 
As shown in Fig. 6, there is almost no charge sitting on 
K sites. This resultpis consistent with previous calcula- 
tion on Na2BagGe4 6 B, which found a rather simple charge 
transfer from Na to the Si skeleton. In their calculation, 
some hybridization are found between Bp. and Si since 
Ba atoms has some low lying 5d orbitalsO. This differ- 
ence also corresponds to the DOS and other electronic 
properties. A high DOS peak at Fermi energy is found 
for Na2Ba 6 Ge460 and this material is superconducting)! 
In comparison, the DOS at Fermi level for KsGe46 is 
moderate and it is not superconducting. 

V. CONCLUSIONS 

We have used first principles SCF pseudopotential 
method to investigate the structural and electronic prop- 
erties of Ge46 and KgGe46 clathrates. The main conclu- 
sion of this work can be made as follows: 

(1) Germanium clathrate Ge46 is found to be a locally 
stable structure with its energy only slight higher than 
that of diamond phase, and its atomic volume is about 
13% larger than diamond phase. 

(2) Ge46 clathrate shows an indirect band gap along 
T-X direction that is about 1 eV higher than the band 
gap of diamond phase. The pentagonal rings in the 
clathrate structure cause the valence band structure of 
Ge46 clathrate to be similar to that of Ge24 fullerene 
cage. The open covalent network structures contribute 
to the large band gaps. 

(3) The K§Ge46 clathrate is metallic with a moderate 
density of states. The valence band structures and DOS 
are similar to those of the pure Ge46, while the conduc- 
tion bands are modified due to the K dopants. Almost 
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